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1. Introduction 

Many phenomena of electrical conduction in small mesoscopic conductors 
have been successfully explained within the framework of the scattering 
approach [1]. The main emphasis of this work is an extension of this ap- 
proach to time-dependent phenomena. Of particular interest are the basic 
requirements which a dynamic conductance of a mesoscopic conductor has 
to satisfy. There is a gap in the way dynamic conduction is treated in a large 
body of the physics literature and in more applied discussions based on sim- 
ple circuit theory. For an electric network we know that its time-dependent 
behavior is crucially determined by its capacitive and inductive elements. 
In contrast to this, time-dependent phenomena are often treated by a re- 
sponse theory of non-interacting carriers. Unfortunately, this approach to 
ac-conductance which has been employed by some of the leading condensed 
matter physicists [2] is prevalent. We emphasize that in order to obtain a 
reasonable answer it is not sufficient to consider non-interacting electrons, 
a Fermi liquid or even a Luttinger liquid with short range interactions. An 
analysis based on an electric circuit model gives an answer which conserves 
the total charge and which has the property that all frequency dependent 
currents at the input and output nodes of such a circuit add up to zero. 
To obtain results which are charge and current conserving, it is necessary 
to consider the implication of the long range Coulomb interaction. This 
leads to a theory which can be considered to be an extension of the work 
of Pines and Nozieres [3] and Martin [4] on bulk systems to mesoscopic 
conductors. Below we present a discussion [5] which for simplicity focuses 
on the low-frequency transport. The benefit of this restriction is its consid- 
erable generality, and that it can be used to analyze the self-consistency of 
the dc scattering-approach (Sect. 3), the leading order nonlinearities of the 
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Figure 1. Mesoscopic capacitor connected to contacts with electrochemical potentials 
/Lti and //2 and electrostatic potentials J7i and U2- After Ref. 7. 

current voltage characteristic (Sect. 4), and the low-frequency dynamical 
conductance (Sect. 5). 

For a two-terminal conductor the low-frequency admittance is of the 
form 



where G(0) is the dc-conductance and -E^ is called the emittance [5]. The 
emittance describes the current response in the leading order with respect 
to frequency and can be associated with the displacement charge passing 
a contact. For conductors with poor transmission as, e.g., a condenser or 
a metallic diffusive wire, the omittance is positive: the current follows the 
voltage as is characteristic of a capacitive response. On the other hand, it 
turns out that the emittance is negative in samples with good transmis- 
sion. The current leads the voltage as is characteristic of circuits with an 
inductance. 

To illustrate the physics contained in this response coefficient E^, we 
present now a few results. Their detailed derivation is presented in Sect. 6. 

Consider first the case of a mesoscopic condenser [6, 7] as shown in 
Fig. 1 with a geometrical capacitance Cq. The measured electrochemical 
capacitance is determined by the geometrical capacitance in series with 
quantum capacitances e^dNi/dE, e'^dN2/dE. Here dNi/dE and dN2/dE 
are the total densities of states of the mesoscopic plates in the regions into 
which the electrical field penetrates. Thus the electrochemical capacitance 
of a mesoscopic condenser is [7] 



G{u;) = G(0) - iiuEf, + 0(a;^) . 



(1) 



C-^ = + (e^dNi/dE)-^ + {e'^dN2/dE)-\ 



(2) 
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Figure 2. Mesoscopic metallic conductor connecting two reservoirs. The low-frequency 
response is capacitive and proportional to 1/6 of its total density of states. 

As it must be, in the macroscopic limit, where the densities of states di- 
verge, one finds Ef^ = Cq. On the other hand, if the density of states at 
the surfaces of the plates are sufficiently small, the geometrical capaci- 
tance can be neglected. For the symmetric case this leads to the emittance 

= (l/4)e'^{dN/dE), where dN/dE = dNi/dE + dN2/dE is the total 
density of states of the two surfaces. 

Next consider a sample which permits transmission of carriers between 
the contacts. To be definite, first consider a metallic diffusive wire con- 
necting two reservoirs (see Fig. 2). Experiments have been carried out by 
Pieper and Price [8] and related theoretical work can be found in Refs. 
[9, 10, 11, 12]. We determine the local potential with the help of Thomas- 
Fermi screening (a local charge-neutrality condition). In the presence of 
an applied dc-voltage, the ensemble averaged potential drops linearly. The 
total density of states in the volume over which the electric field is non- 
vanishing is denoted by dN/dE. We find that the emittance of such a 
metallic diff'usive wire is positive like a capacitance and is given by [13] 

E^ = {l/Q)e\dN/dE) . (3) 

It is not the total density of states that counts! Without screening and in 
the presence of a dc potential-drop, one-half of the states would be filled. 
Screening, however, reduces the naively expected result (l/2){dN/dE) by 
a factor 1/3. Such a 1/3-reduction is familiar from the theory of shot noise 
of metallic diffusive conductors were it is found that the actual shot noise 
is 1/3 of the full Poisson noise [14, 15, 16]. The physical origin of the factor 
1/3 in these two problems is quite different but from a purely mathematical 
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point of view the origin [15] of the factor 1/3 is the same. 

Next consider a balhstic wire of length L. Now, for dc transport the 
ensemble-averaged potential drops as wc enter the wire from the reservoir, 
it is constant in the wire, and it drops again as we leave the wire. We find 
then a negative (inductive) emittance 



where dN/ dE is the total density of states of the wire in the volume where 
the potential is uniform. 

As a further example, consider a resonant double barrier [17, 18, 19]. 
Suppose for simplicity that it is symmetric and suppose that it is reasonable 
to determine the potential in the well from a charge neutrality condition. 
If we denote by T, i? = 1 — T, and dN/dE the transmission probability, 
the reflection probability, and the density of states in the well, respectively, 
the emittance can be written in the form [19, 20] 



Note that Eq. (5) interpolates between the emittance of a symmetric con- 
denser with large geometrical capacitance and a ballistic wire. The emit- 
tance is negative (inductive) at resonance T = 1 and crosses zero at T = 
R = \/2 where the Fermi level is a half- width of the resonant level above 
(or below) the energy of the resonant state. If the Fermi level is further 
than a half-width away from the resonant level the emittance is positive 
(capacitive). The case of a double barrier which is asymmetric has been 
considered in Ref. [19]. A result for single tunneling barriers which goes 
beyond Thomas- Fermi screening will be discussed elsewhere [20] . 

Finally, consider a two-dimensional electron gas patterned into the 
shape of a Corbino disk (Fig. 3a) or a Hall bar (Fig. 3b) with two metallic 
contacts. In the range of magnetic fields over which the Hall conductance 
is quantized the dc-resistances can be evaluated solely by considering edge 
states [21]. In a non-interacting theory the edge states intercept the Fermi 
energy with a finite slope and have a local density of states per unit length 
given by dn{s)/dE = l/hv{s), where s is a coordinate along the edge state 
and where v{s) is the velocity of carriers at point s. The integrated density 
of states along the edge channel k is thus dN^/dE = f ds{dnk/dE). Sup- 
pose for simplicity that the magnetic field is such that we have only one 
pair of edge states. The Corbino disk acts then as an insulator with vanish- 
ing dc conductance G(0) = 0. The long-range Coulomb interaction between 
the inner and the outer edge state can be described by a geometrical ca- 
pacitance Co (which depends logarithmically on the width of the sample) . 
As one expects, we find that the Corbino-disk exhibits an emittance which 



Ei, = -{l/4)e\dN/dE) 



(4) 



= (l/4)(i? - T)e'^ (dN/dE) . 



(5) 
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Figure 3. Hall conductors: (a) Corbino disk, (b) Hall bar. After Ref. 22. 



equals the cleetrochemical capaeitanee [22] iJ^j = C^^, where is given by 
Eq. (2). A topologically different conductor of the same width and length 
but with its contacts arranged at its ends is shown in Fig. 3b. Now the 
two-terminal dc-conductance is quantized and given by G(0) = /h for 
spin-split Landau levels. Obviously, the low frequency response of this con- 
ductor is not capacitive but dominated by transmission of carriers between 
different reservoirs. In particular, we find the inductive-like emittance [22] 
= —Cfj,. Thus the ac-response of a quantized Hall conductor is de- 
termined by the way the edge states are connected to reservoirs. In the 
Corbino disk each edge state returns to the reservoir from which it em- 
anates. In contrast, in the two-terminal Hall-bar, each edge state connects 
to a reservoir that is different from the one from which it emanates. We do 
not discuss this example anymore in the sequel. A forthcoming publication 
[22] presents a general formula for the emittance which is applicable to a 
wide variety of samples with different edge-state topologies. 

2. The Scattering Approach 

Consider a conductor connected to a number of contacts [23]. The contacts 
are labeled with the greek indices a. We assume that the distance between 
these contacts is so small that transmission from one contact to another one 
can be considered to be phase coherent. Thus we assume that scattering 
inside the conductor is purely elastic. For small deviations of the electro- 
chemical potentials away from their equilibrium value the dc-current la at 
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probe a is given by [23] 




Here, is the number of quantum channels with thresholds below the 
equilibrium electrochemical potential in contact a, and / is the Fermi func- 
tion. Carriers in the A^^ incident channels have a total combined probability 
Raa for reflection back into contact a. Carriers incident in contact /3 have 
a total probability Tq,^ to traverse the sample into contact a. Equation 
(6) is a quantum-mechanical Kirchhoff law. It states the conservation of 
the current at an arbitrary intersection of a mesoscopic wire and that the 
currents are only a function of voltage differences. In the present context 
these features are a consequence of the unitarity of the scattering matrix S 
and its behavior under time reversal. The fact that the scattering matrix 
is unitary and that the microscopic equations are reversible implies that 
under a reversal of magnetic field the scattering matrix has the symme- 
try S^(-B) = S(— S). The scattering matrix S for the conductor can be 
arranged such that it is composed of sub-matrices Sq,^ which relate the in- 
cident current amplitudes in probe /3 to the outgoing current amplitudes in 
probe a. In terms of these scattering matrices the probabilities introduced 
above are simply T^p = Tr(s^^SQ,^) if a 7^ /3 and Raa = 7^'^(sJ,qSq,q). The 
trace is over the quantum channels in each probe. As a consequence of the 
unitary properties of the scattering matrix, we have = Raa + Jl/sTafS 
and Na = Raa + J^p Tf3a which imply current conservation J^a = 0. Fur- 
thermore, the currents depend only on differences of the electrochemical 
potentials. 

Let us restate these properties in terms of conductance coefficients 
Gap = la/Vjs. Taking into account that in a reservoir the voltages and 
electrochemical potentials move in synchronism, = eV^, we have [23] 

G„„{0) = ^/<iE(-£){W.-ii„.) = £!/<i£(-£)i;T„, . (7) 

for q; = /3 and 

Gam = ~JdEi-^)T^p . (8) 
Current conservation means for the conductance matrix 



a 



(9) 
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for any /?. On the other hand, the statement that the currents can only 
depend on voltage differences leads to 



We emphasize the simple properties (9) and (10) since later on we point 
out that they are also valid for Ga/3{(^)- 

Before proceeding wc briefly comment on closely related work. In his 
work Landauer determined the voltage drop across an obstacle with the 
help of a local charge neutrality argument applied to the perfect sections of 
the conductor on either side of an obstacle [24] . For an obstacle with trans- 
mission probability T, this leads to a resistance proportional to (1 — T)/T. 
This resistance has the obviously correct property that it vanishes in the 
limit of perfect transmission. Already Engquist and Anderson pointed out 
that voltages measured at contacts are determined not by a charge neutral- 
ity condition but by adjusting the electrochemical potential of a voltage 
probe in such a manner that the voltmeter draws zero net current [25]. 
Due to the special geometry considered, the restriction to a one channel 
conductor, and the restriction to a phase-incoherent voltage measurement 
[26] , the answer found by Engquist and Anderson is the same as that found 
by Landauer. The key notion that the resistance across an obstacle should 
vanish if the obstacle permits perfect transmission persisted long after the 
work of Engquist and Anderson. The task seemed, therefore, to be to find 
resistance formulae [27] of the type (1 — T)/T. In a work that largely centers 
on this notion, Imry [28] observed that if one considers in the discussion 
of Engquist and Anderson not the voltage drop across the obstacle but the 
voltage drop between the current source and sink one obtains a resistance 
that is proportional to T~^, i. e. a conductance that is proportional to T. 
Equation (6), in contrast, was motivated by experiments by Benoit et al. 
[29] in which the role of the current and voltage source was exchanged. 
Such experiments imply that we need a formulation of electrical transport 
in which apriori all contacts are treated equivalently and on an equal footing 
[23]. In a specific arrangement of current and voltage sources the voltages 
are a posteriori determined by a zero current condition. This leads to con- 
ductances which are electrochemical quantities and which obey the basic 
requirements of a transport theory. These basic requirements include a reci- 
procity symmetry [23, 29, 30] (which is a consequence of microreversibility 
and the irreversibility of thermodynamic electron reservoirs) and include a 
fluctuation-dissipation theorem [31]. Below we discuss in more detail the 
electrochemical requirements which lead to Eq. (6). 

It is sometimes implied in the theoretical literature, that the derivation 
of conductance formulae is just a question of applying formal linear response 




(10) 
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Figure 4- Electrochemical and electric potentials of a mesoscopic capacitor: An in- 
crease of the electrochemical potential /Lti by d//i causes the band bottom to shift by 
edU{r) = ui{r)dni where wi is the characteristic potential of contact 1. After Ref. 5. 

theory correctly to this problem. This is not borne out by the history of this 
field: Depending on one's preconceived notions quite different results can 
be derived [32]. But it is certainly true, that the derivation of these results 
from linear response has made them more acceptable. This is, however, 
neither a question of rigour nor depth but simply a consequence of the fact 
that even simple results arc often only accepted if they arc embellished 
by a sufficiently complicated derivation. For very concise linear response 
discussions of Eq. (6) and some still open related questions we refer the 
reader to Ref. [33, 34, 35]. 

3. Potential Distribution in a Mesoscopic Conductor 

The scattering approach as discussed above seems like a simple black-box 
approach. If we know through which contact carriers leave the sample and 
know the current that is injected by a reservoir we can find the total cur- 
rents. The issue which is not trivial is the fact that the voltages at the 
contacts must be well-defined in order for Eq. (7) and Eq. (8) to be the 
correct answer. In many articles on the subject one finds a conductance 
attributed to a wire which is strictly one-dimensional, or one finds pictures 
of mesoscopic conductors with leads that are narrower then the mesoscopic 
sample itself and are called reservoirs. Such geometries, as we now show, 
do in fact not lead to the conductances given above. This criticism is not 
novel: it has been aired in a number of papers by Landauer [36] and one 
of the authors [31] - unfortunately without much success. The following 
discussion is most closely related to work by Levinson [37] but in detail 
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follows Ref. [5]. 

For Eqs. (7) and (8) to be valid, the electric potential Ua in the reservoirs 

a must follow the electrochemical potential fi^ in this contact. To generate 
transport we must consider a non-equilibrium situation. The electrochemi- 
cal potentials in the contacts must be allowed to differ from the equilib- 
rium chemical potential /iq. Suppose that the increment of the electrochem- 
ical potential in contact a is 5fia = fJ-a— /^o- The electric potential changes 
from its equilibrium configuration Ueq{r) to a new non-equilibrium configu- 
ration i7([jLta], r). Here the argument [^a\ indicates that the non-equilibrium 
potential depends on the electrochemical potentials in the contacts. Equa- 
tions (7) and (8) presuppose that deep in contact a the electrochemical 
and electric potential change in synchronism, i. e. Sjia = e^U{[iJ,a],r). This 
is a consequence of the charge neutrality deep in the reservoirs. The differ- 
ence between the electrochemical potential and the electrostatic potential 
is the Fermi energy, Epaij^) = l^a — eU{r). This chemical potential Ep^ 
determines the charge in the neighborhood of r. Deep in the reservoir the 
charge cannot change, even if we bias the conductor. Consequently we have 
= SEpai^^) = SfJ'a — G^U{[na], r). Notc that wc have taken here the chem- 
ical potential to be space depended. However, if the quantum coherence 
of the wave functions is taken serious then also deep in the reservoir, i. 
e. in the wide wire, a single energy Epa is all that is needed to specify 
the chemical potential. The energy dispersions in a wide wire are of the 
type Ean{k) = h'^k'^/2m + E^.^^ + eUa- Here h^k'^/2m is the longitudinal 
energy for motion along the wire, i?^^ is the energy for transverse motion 
(the channel threshold) and eUa is the equilibrium potential (the bottom 
of the conduction band) in contact a. The energy dispersion depends in 
an explicit manner only on a single spatially independent constant eUa- 
Instead of the spatially dependent relationship between the electrochemi- 
cal, chemical, and electrostatic potentials we have for coherent conductors 
the relation Epa = l^-a — eUa- In either case the differential relationship is 
spatially independent in the lead. In order to discuss the validity of Eqs. 
(6)- (8), we have thus to find the electrostatic potential and the conditions 
under which it changes synchronously with the electrochemical potential in 
the reservoirs. 



3.1. CHARACTERISTIC POTENTIALS 

The electrostatic potential J7([/Lia],r) for mesoscopic conductors is a func- 
tion of the electrochemical potentials of the contacts, and a complicated 
function of position. Small increases in the electrochemical potentials 
bring the conductor to a new state (see Fig. 4) with an electrostatic poten- 
tial U{[iJLa + 5A*a],r). The difference 5U between these two potentials can 
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be expanded in powers of the increment in the electrochemical potential. 
To linear order we have 

e6U{[na],r) ='^Ua{r) Sfia ■ (H) 

a 

Here, Ua{r) = edU{[fj,a],i'')/dfia\fMa=f^o^ with a = 1,2 are the characteristic 
potentials [5] which determine the electrostatic potential inside the sample 
in response to a variation of an electrochemical potential at a contact. 

Suppose for a moment that we increase all electrochemical potentials 
simultaneously and by the same amount, Sfi^ = Both before and after 
the change the conductor is at equilibrium, hence the physical properties 
of the conductor remain unchanged. 

Consequently, the shift of the electrochemical potentials must be ac- 
companied by a shift e6U = 6fi of the electrical potential. Imposing this 
condition on Eq. (11) implies that the sum of all characteristic potentials 
is equal to one at every space point [5], 

J2ua{r)^l. (12) 

a 

Equation (12) is a consequence of the long-range Coulomb interaction. It 

is the most important result of this work: the conservation of charge under 
the application of a dc or ac bias and the conservation of current are, as we 
will show, a consequence of Eq. (12). 

Let us now return to our original problem and consider what happens 
if we increase just one electrochemical potential, say in reservoir a, by 
6iJ,a- Obviously, the condition that the electrochemical potential and the 
electrostatic potential move in synchronism deep inside reservoir a implies 
that the characteristic function Ua{r) = 1 for r deep inside reservoir a. 
Together with Eq. (11), this implies that Eqs. (6)-(8) are valid if and only 
if the characteristic potentials have the property that tta(r) = 1 for r deep 
in contact a and no,(r) = for r deep in any other contact. 

The electrostatic potentials are determined by the charge distribution 
in the sample. As we increase the chemical potential of contact a keeping 
all electrostatic potentials fixed, the additional charge 

Sn{r) = {dn{r , a) / (IE) 6 (13) 

enters the conductor. Here, dn{r,a)/dE is the injectivity of contact a into 
point r of the sample. With the help of the scattering states *an(r) which 
have unit amplitude in the incident channel n in lead a, the injectivity can 
be expressed as 

dn{r,a)/dE = ^(/ii;an)-'|*an(r)P (14) 

n 



ELECTRICAL CONDUCTION 



11 



where Van is the velocity of carriers at the Fermi energy in channel n in 
contact a. Equation (13) gives of course not the true density variation. 
The injected charges induce a change in the electrostatic potential which 
in turn implies an induced contribution to the density, Sriind, which has to 
be determined self-consistently. The total charge density is 

6n{r) = {dn(y,a)/dE) Sna. + (^nmd(r) ■ (15) 

The induced charge density is connected to the electrostatic potential via 
the response function n(r,r') (Lindhard function): 

8nind{v) = - j d\' n(r,r') e5U{r') . (16) 

The response function can be expressed in terms of the scattering states 
(Green's function of the Schrodinger equation). Note that the Lindhard 
function describes the variation of the charge density not only of the mobile 
electrons which can be reached from the contacts but also of the localized 
states. For the purpose of our discussion we simply assume that this re- 
sponse function has been calculated and is known. There is one property 
of the Lindhard function which is needed later on and which is a simple 
consequence of the invariance of the electrical system under a global poten- 
tial shift. A simultaneous change in all electrochemical potentials injects a 
charge 

5n{v) = ^idn{r, a)/dE)6fi + SrHndir) = . (17) 

a 

Taking into account Eq. (16) and e6U = S/j,, one concludes [5] 

dn{r)/dE = j d^r' n(r,r') . (18) 

The quantity on the left hand side side of Eq. (18), 

dn{r)/dE = Y,{dn{r,a)ldE) , (19) 

a 

is the local density of states which equals the sum of all inject ivities. Equa- 
tion (18) connects a chemical response quantity, the local density of states, 
to the Lindhard function 11. It is in this regard similar to the Einstein re- 
lation between a diffusion constant and a conductivity. 

Now we come back to the case where the voltage is changed only 
in contact a. By inserting Eq. (15) into Poisson's equation and taking 
e5C/(r) = Ua{r)5iJLa into account, we find that the characteristic potential 
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lia(r) is the solution of a field equation with a non-local screening kernel 
and a source term given by the injectivity of contact a, 

-Aua{r) +4Tre'^ j d^r'n{r,r')ua{r') = Ane^{dn{r,a)/dE) . (20) 

We define the Green's function ^(r,ro) as the solution of Eq. (20) with the 

source term e dn{r, a)/dE replaced by a localized test charge eS{r — rg) at 
point tq. The characteristic potential ita(r) can then be written in the form 

uair) = J d\'g{v,v') {dn{v',a)/dE) . (21) 

Using Eqs. (12) and (19), a summation over a implies for Green's function 
the property [5] 

/ d\'g{y,v')Y,{dn{v',a)/dE) = [ d\'g{y,v') {dn{v')/dE) = 1 . (22) 

The same relationship follows from the condition that the sum of all in- 
duced charge densities plus the test charge is zero. 

Now we find the condition for the electrical self-consistency of Eqs (6). 
According to Eq. (22) the characteristic potential is equal to unity if the 
Green's function is convoluted with the local density of states. Therefore, 
we must have that the injectivity dn{r,a)/dE deep in contact a is equal 
to the local density of states dn{r)/dE. This requires that nearly all (in 
a thermodynamic sense) electrons approaching the contact a be reflected 
into the reservoir. If the conductor and the reservoir consist of the same 
material then the reservoir must be wide compared to the mesoscopic con- 
ductor. In semiconductor samples with metallic contacts, on the other hand, 
the contact might be actually narrow compared to the dimensions of the 
semiconductor since the density of states at the Fermi energy of the metal 
is much larger then that of the semiconductor. This is the case, for in- 
stance, in Ga/As-samplcs used to measure the quantized Hall effect [21]. 
Our emphasis that self-consistency requires geometries which are of a wide- 
narrow-wide geometry deserves further discussion: the notion that a por- 
tion of length L of a purely one-dimensional conductor has a conductance 
G = {e^ /h)T (per spin) seems to be widely accepted. In contrast, from the 
point of view taken here a strictly one-dimensional conductor cannot be 
characterized by a conductance. 

3.2. CHARGE CONSERVATION 

Let us now demonstrate that if the conditions of self-consistency are met, 
then the application of an external bias to the conductor preserves the 
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total charge in the system. To be more precise, imagine a volume which 
encloses the entire conductor including a portion of the reservoirs which 
is so large that at the place were the surface of Q intersects the reservoir 
all the characteristic potentials are either zero or unity. We demonstrate 
charge conservation for the case that we increase jUi by a small value S/ii 
above the equilibrium chemical potential. The change in density is given 
by Eq. (15). Expressing the induced charge with the help of the Lindhard 
function n(r, r') and using the Green's function to relate the characteristic 
potentials to the injected density gives 

^ - I MVU(r. r'W. r")^^) *m (23) 

The total variation in charge is Jq d^r dn{r) = 

where we used Eq. (18) and the symmetry property n(r, r') = n(r',r) of 

the Lindhard function. Equation (22) then implies that the total variation 
in charge inside the volume Q vanishes. This can be understood in the 
following way. According to the law of Gauss the charge included in a vol- 
ume O is Jqq EriS = 4:TtQ. The charge in CI is conserved if the electric flux 
through the surface of vanishes. This means that for this conductor all 
electric field-lines which are generated when we bias the sample have their 
sources and sinks within the volume CI. Application of a bias voltage to 
an electrical conductor results in a redistribution of the charge within our 
sample but not in an overall change of the charge. If the conductor is poor, 
i. e. nearly an insulator, the reservoirs act like plates of capacitors. In this 
case long-range fields exist which run from one reservoir to the other and 
from a reservoir to a portion of the conductor. But if we chose the volume 
O to be large enough then all field lines stay within this volume. 

In all these considerations we have implicitly assumed that our conduc- 
tor and the reservoirs is all that counts. Such a situation might be realized 
for the metallic mesoscopic structures fabricated on insulating substrates. 
But this picture is certainly not complete if we deal with modern meso- 
scopic semiconductor structures which are often defined with the help of a 
number of nearby gates. In such a case we must take a broader view and 
include inside our volume Cl not only the conductor of interest but also 
all of the nearby gates [5]. ^Prom an electrostatic point of view, we deal 
then not only with the mesoscopic object of interest but we have to take 
into account the nearby electrical bodies used to define this object. In such 
a case the overall charge is still conserved, even though the total charge 
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on the mesoscopic conductor of interest varies with the apphed bias. The 
theory presented here can easily be extended to this case [5] . 

4. Nonlinear I-V Characteristic 

As an apphcation of the discussion given above, let us consider the non- 
hnear 1-V characteristic. The discussion presented here can also be carried 
out in terms of an external and internal response and provides a nice illus- 
tration of these concepts [5] . Nonlinearities in metallic mesoscopic samples 
have been analyzed by Al'tshuler and Khmelnitskii [38] using diagrammatic 
techniques without a self-consistent potential. For transmission through a 
tunnel contact the effect of a potential which changes with increasing ap- 
plied bias has been investigated by Prenkel [39]. Landauer has pointed to 
the necessity of a self-consistent treatment of the internal potential [40]. 
In this section we derive the current-voltage characteristic taking into ac- 
count nonlinearities which are a consequence of the increase of the exter- 
nal electrochemical potential differences as well as the changing internal 
(self-consistent) potential distribution. We focus on the leading nonlinear 
correction of the low-voltage ohmic behavior of the sample. The I-V charac- 
teristic of a mesoscopic sample is in general rectifying, i. e. I{V) ^ —I{—V). 
Furthermore, since rectification also depends on the internal potential and 
since the internal potential in conductor k depends on the charge distri- 
bution of other nearby conductors, the rectification properties of a small 
sample are dependent on its entire electric environment. Nevertheless, we 
consider for simplicity a conductor which is in electric isolation. Reference 
[5] presents a more general result being valid if there are additional nearby 
conductors like gates or capacitors. 

To proceed wc view the scattering matrices as a functional of the po- 
tential distribution Sai3{E,U{[iJ,a],T^)) and expand U away from the equi- 
librium potential-distribution. The scattering matrix in the neighborhood 
of the equilibrium reference-state (index 0) is Sa/3{E,eU{r)) = s°^(£') -|- 
/ d'^r'{6s^p/e5U{r'))e5U{r'). Here, eSU{r') can be expressed in terms of 
the characteristic potentials and the electrochemical potentials of the reser- 
voirs. The total current at probe a can be found by the same considerations 
that lead to Eq. (6). However, we stop short of linearizing the resulting 
expression in the electrochemical potentials. The neglect of any inelastic 
scattering in the presence of large applied voltages is of course a limitation, 
but since we only focus on the quadratic term in the voltages this limitation 
might be not be so serious. The current is [31] 

= (e/h) ^ J dEfpTr [l^^,^ - slp{E, U{v))s^^{E, U{v))\ , (25) 
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where fp is the Fermi function belonging to reservoir /?. The sum over 
all currents at all terminals is still zero due to the unitarity of the scat- 
tering matrix, hence the current is conserved. In order that the current 
depends on voltage differences only, it is necessary to treat the potentials 
self-consistently. We expand Eq. (25) in powers of the electrochemical po- 
tential deviations Sfia = eV^, 



^a = Yl dafiVp + (1/2) gap-yVpV^ 

p Pi 



(26) 



The terms linear in the electrochemical potentials are determined by the dc- 
conductances ga/3 = (e^/h) J {—df / dE)Tr lap — Sa^(-E')sa/3(-E)j which are 
a functional of the equilibrium reference-potential only. The leading order 
nonlinear terms are given by transport coefficients which are composed of 
an external and an internal response: 



(27) 



The external response arises from the expansion of the Fermi functions in 
powers of the electrochemical potentials and is given by [5] 



Sa/S-y 



e 

T 



J dE{- 



dh 
dE 



)Tr 



dE 



+ 



t 

a/3 



dE 



Sa/3 



(28) 



where /o is the equilibrium Fermi function. Let us examine the external 
response for the case that we have a two-terminal conductor. ^Prom the 
unitarity of the scattering matrix we find that all the non- vanishing second- 
order conductance coefficients are equal in magnitude, 



6 6 6 

9lll = ~9l22 = "5211 



5222 



(29) 



Furthermore, these coefficients can also be expressed just as the energy 
derivative of the transmission probability, 



9in = jJdEi--)- 



(30) 



Thus up to the second order the contribution of the external response to 
the current is 



-h = giVi - V2) + {l/2)glr,{V,^ Vi) 



(31) 



Despite the fact that currents are conserved this is an unphysical result. 
The quadratic term depends not only on the voltage difference but on the 
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individual voltages. Equation (31) would predict that we should observe a 
different current depending on whether we rise the voltage of the left con- 
tact by Sfii = eV\ or whether we decrease the voltage on the right contact 
by = — eVi as compared to the equilibrium value of the electrochemical 
potential. This is a simple example which demonstrates why the calcula- 
tion of a nonlinear current voltage characteristic without the self-consistent 
adjustment of the electrostatic potential makes no sense. 

Let us now consider the internal response. The internal response is a 
consequence of the change in the potential distribution and is given by [5] 



Uf,{r). (32) 



Note that the internal response contributes only to quadratic order in the 
voltage. The linear conductance is a purely external response. If we now 
add external and internal response, take into account Eq. (12) and that the 

integral over Q of an internal response term with n = 1 is equal to minus 
the external response with the functional derivative replaced by an energy 
derivative, we find [5] 

fl'a/37 = gafit = 13 Safi'y = . (33) 
a p J 

For a two-terminal conductor the second-order conductance coefficients 
obey 5111 = — 5112 = -5121 = 5122 = —5211 = 5212 = —5221 = —5222 
Consequently, the currents are 

h = -h = 9{Vi - V2) + (l/2)5iii(Fi - V2f + ... . (34) 

Now, the current depends only on the voltage difference as it must be. In 
contrast to the external response which could simply be expressed in terms 
of energy derivatives of the transmission probability, the total response de- 
pends on the charge distribution inside the conductor. 

As a simple example, we consider an asymmetric resonant double bar- 
rier. The long lived state has a decay width Fi to the left and F2 to the 
right. For simplicity, assume that the potential in the well is determined 
by a local charge-neutrality argTimcnt (sec Sec. 6). The characteristic po- 
tentials in the well are ui = Fi/F and U2 = r2/r, where F = Pi -|- r2 is 
the decay width. ^From Eqs. (28) and (32) we find for the second order 
conductance coefficient giw = {e^ / h){dT /dE){l — 2ui) and hence [13, 20] 

ffiii = ieyh){dT/dE){r2 - Fi)/F . (35) 
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In summary, we emphasize that the nonhnearity cannot be discussed with- 
out a concern for the way the potential drops in the interior of the conduc- 
tor. 

5. Frequency Dependent Conductance 
5.1. EXTERNAL RESPONSE 

Wc arc interested in the dynamical response of the conductor. A time- 
dependent voltage oc exp(— iwt) can be applied across two terminals, be- 
tween a terminal and a nearby gate, or between two nearby gates. We want 
to know the currents which appear as a consequence of these oscillating 
voltages at the contacts of the conductor or at the contacts to the nearby 
gates. We are seeking the admittance matrix 

Gapiuj) = Ia{^)/Vp{u) . (36) 

Again we consider for simplicity the case of a two-terminal conductor in 
electrical isolation. The case where there are a number of nearby gates 
(capacitors) or other conductors has been the subject of a number of dis- 
cussions [5, 46]. Thus the indices a, j3 take the values 1 and 2 for the left 
and right contact, respectively. If there are no other nearby electrical con- 
ductors then all electric field lines which emanate from the conductor also 
return to the conductor or to the reservoirs. Even in the presence of time- 
dependent voltages applied to this conductor the reservoirs remain locally 
charge neutral, i. e. the electric field lines emanate from a reservoir only in 
the region where the transition to the conductor occurs. If the conductor 
is very short then the reservoirs act like the plates of a capacitor and, due 
to long range Coulomb intercation, field lines connect the surfaces of the 
two reservoirs facing each other. As in the dc-case there exists, therefore, a 
volume n which is so large that there is no electric flux through its surface. 
Consequently, if we include all components of the system within Q, then the 
total charge within this volume is zero, i. e. all currents at the terminals 
must add up to zero. Furthermore, since a potential which is uniform over 
the entire volume Q, is of no physical consequence the resulting currents 
must depend on the potential differences only. Therefore, Eqs. (9) and (10) 
hold also for the dynamic conductance: the rows and columns of the dy- 
namic conductance matrix Gapi^) must add up to zero. We call such a 
discussion of the ac-transport a charge and current conserving theory. Be- 
low we illustrate the features of such a theory for the case of low frequencies 
only. But the extension to a larger range of frequencies and to nonlinearities 
must follow the very same line of thought. 
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Figure 5. Decomposition of the total density of states of a two terminal conductor. 
After Ref. 13. 



5.2. DECOMPOSITION OF THE DENSITY OF STATES 

The total density of states in the conductor inside the volume is a sum 
of four contributions [34] , 



dN/dE = Y,{dNa[5/dE) , 

a/3 



(37) 



where 



dN, 



af3 



dE 



Am 



Tr 



ds, 



•al3 



•^1^ dE 



dE 



(38) 



are the partial densities of states. Fig. 5 gives a schematic representation of 
the partial density of states. The partial density of states dNu/dE consists 
of carriers that originate in contact 1 and return to contact 1. The partial 
density of states dN2i/dE consists of carriers that originate in contact 1 
and are transmitted to contact 2. It turns out that the external response 
is determined exactly by these four partial densities of states. The partial 
densities of states represent a decomposition of the total density of states 
both with respect to the origin of the carriers (injecting contact, right in- 
dex) and the final destination of the carriers (emitting contact, left index). 

At constant electrostatic potential the total charge injected into the 
conductor under a simultaneous and equal increase of the chemical poten- 
tials at its contacts SEpi = SEf2 = eSji is 5Q^ = eY^^p{dNai3/dE)5iJ,. 
To find the current at contact a we need to know which portion of this 
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charge enters or leaves the conductor through contact a, i. e. how the to- 
tal charge is partitioned on the two contacts. The answer to this question 
was found by one of the authors, Pretre and Thomas [34] using a linear 
response calculation. The following simple argument leads to the same re- 
sult. The scattering matrix determines the current amplitudes of the 
outgoing waves in contact a as a function of the current ampUtudcs of the 
incident waves in contact /3. The charge which is injected by an increase 
of the Fermi energy at contact 1 is (edNu/dE + edN2i/dE)eVi{uj). Only 
the additional charge 5Qi{w) = e{dNn/dE)eVi{uj) leads to a current at 
contact 1, whereas SQ2{uj) = e{dN 21 / dE)eVi{uj) is determined by carriers 
which leave the conductor through contact 2. Therefore, a variation of the 
Fermi levels of the contacts (3 causes the current 

6Il{uj) = -iue^ Y.{dN^p/dE)Vf^{u;) (39) 
P 

at contact a. Since direct transmission between contact 1 and 2 is possible, 
an oscillating voltage causes in addition at these contacts a current de- 
termined by the dc-conductance. Thus the leading external low-frequency 
current response to an oscillating chemical potential eVa{uj) is given by [34] 

Gl^{Lu)=Gap{0)-iue\dN^p/dE) . (40) 

This external response is not current conserving. Since the dc-conductan- 
ces satisfy J2aGal3{^) = Jl/s Ga^iO) = 0, one finds that to leading order 
in frequency Z)qG'q^(w) is proportional to the total charge injected from 
contact P into the conductor. The injected charges create an internal, time- 
dependent electric potential SU{r, t) which in turn causes additional cur- 
rents. In the following subsection, we investigate the response to such an 
internal electrostatic potential. It turns out that this requires a detailed 
knowledge of the charge distribution. 

The density of states (37) has been obtained assuming a perturbation 
which is (in a mathematical sense) asymptotically far away from the sam- 
ple. The derivative with respect to the energy is a consequence of the 
asymptotic nature of this perturbation. Physically what counts is the den- 
sity of states in a finite volume fi. To obtain these densities it is better to 
first calculate the local densities corresponding to Eq. (38) and to integrate 
these local densities over the volume Q. As shown below the local densities 
are not given by energy derivatives of the scattering matrix but by func- 
tional derivatives with respect to the local potential eU(r). The densities 
determined by integration of such local densities of states differ in general 
from a simple energy derivative by a quantum correction [41, 42]. The dif- 
ference vanishes in the semi-classical (WKB) limit. The same is of course 
valid for the derivatives in Eq. (28). 
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We conclude this section with a remark on locaUzed states. Equation 
(37) is not complete. A conductor might also contain a contribution to the 
density of states from localized states, in addition to the extended scatter- 
ing states considered so far. For the external response the localized states 
play no role, but later when we consider the screening the localized states 
are important. To be brief, however, we do not discuss here their role in 
detail. 

5.3. RESPONSE TO AN OSCILLATING ELECTROSTATIC POTENTIAL 

To investigate the self-consistency of dc-transport we already discussed the 
local charge distribution inside the conductor. In Eq. (13) we introduced 
a local partial density of states which wc called the injectivity (see Fig. 6). 
Now we introduce additional local partial densities of states which permit 
eventually to write the dynamical conductivity in a simple and transparent 
manner. We are interested in the currents generated at the contacts of a 
sample in the presence of an oscillating potential 5U{r,t). We can Fourier 
transform this potential with respect to time and consider a perturbation 
of the form 



Since the potential is real we have U-^ = U^^^. The response to such a po- 
tential can be treated using a scattering approach [19]: due to the oscillating 
internal potential a carrier incident with energy E can gain or loose mod- 
ulation energy hui during reflection at the sample or during transmission 
through the sample. The amplitude of an outgoing wave is a superposition 
of carriers incident at energy E and at the side-band energies, Ezizhcv. In the 
low-frequency limit the amplitudes of the out going waves can be obtained 
by considering the scattering matrix s^piU (r, t),E) to be a slowly varying 
function of the potential U (r, t). Since the deviations of the actual potential 
away from the (time-independent) equilibrium potential Ueq{v) are small, 
we can expand the scattering matrix in powers of 5U (r, t) = U (r, t) — Ueq{r) 
to linear order 



s„^(;7(r, t), E) = s^p{Ueg{r), E) + {5s^p/5U{v))5U{v, t) . (42) 



Here we have introduced the local partial density of states [5, 19] 



5U{r, t) = u(r)(C/_|_tj exp(— iwt) + U-^^ exp{+iu}t)) . 



(41) 



Evaluation of the current at contact a gives [19] 




(43) 



dn{a, r) 



dE 




l^eSUir) eSU{r) 



Sa/3 



(44) 
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Figure 6. Decomposition of the local density of states of a two terminal conductor into 
injectivities. After Ref. 13. 

which we call the emissivity (see Fig. 7). It describes the local density of 
states of carriers at point r which are emitted by the conductor at probe 
a. A more detailed derivation of Eq. (43) can be found in Ref. [19]. It is 
useful to express the response to the internal potential in the from of a 
conductance defined as SP^iuj) = ^^(0))?/+^^. Comparison with Eq. (40) 
gives for the internal conductances 



Below we use this internal response to complete the calculation of the total 
current. Before doing this it seems useful to pause for a moment and to dis- 
cuss in more detail the local density of states which determine the internal 
response. 

5.4. DECOMPOSITION OF THE LOCAL DENSITY OF STATES 

In Eq. (14) we have expressed the injcctivity (13) with the help of scatter- 
ing states. Now we give an expression of the injectivity in terms of deriva- 
tives of the scattering matrix. Expressions which relate wave functions to 
functional derivatives are known from the discussion of the characteristic 
times occurring in tunneling processes [43, 44]. Consider for a moment a 
one-dimensional scattering problem with a potential V{x) in an interval 
(—a, a). The scattering matrices are s^p where a and /3 take the values 
1 and 2 to designate left and right, respectively. Of interest is the time a 
particle dwells in this region irrespective of whether it is ultimately trans- 
mitted or whether it is ultimately reflected. There are two dwell times tdq. 




(45) 
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for the particles arriving from the left or from the right. In terms of the 
scattering states '^a{x) and the incident current I, the dwell time is given 
by [44] 



rx+a 

TDa = / dx\^a{x)\Vl • (46) 
Jx 



To find the time a particle dwells in an interval (,t, ,t + a) an infinitesimal 
uniform perturbation dV is added in this region to the potential V{x). It 
is found that the dwell time is then related to the scattering matrix via the 
following relationship [43, 44] 

TiJa = nimMsia| |-|s2a| 1 • (47) 

where Im denotes the imaginary part. For a plane- wave scattering state 
with wave vector k the current is v = hk/m. Thus a comparison with Eq. 
(14) shows that for a single quantum channel the dwell time is related to 
the injectivity by 

px+a 

TDa/h= / dx (dn{x,a)/dE) . (48) 

Jx 

This means that their exists a simple relationship between local density of 
states and derivatives of the scattering matrix with respect to potentials. It 
is easy to extend this relation to the case of an arbitrary space dependent 
potential and to an arbitrary number of channels. The final result is that 
the injectivity is given by [5, 19] 



dn{r,a)/dE = -^Xl^^ 



^^"e(^C/(r) e5C/(r)^^" 



(49) 



The injectivity contains information about the origin of the particles: it is 
important through which contact the carriers enter the conductor. Conse- 
quently, the summation is over the first index of the scattering matrices. 
In contrast, the emissivity defined in Eq. (44) contains information about 
the future of the carriers: it is important through which contact the carri- 
ers leave the sample. The summation is thus over the second index of the 
scattering matrix. The sum of all the injectivities (see Fig. 6) or the sum 
of all the emissivities (see Fig. 7) is equal to the local density of states. 



dn{r)/dE = Y^dn{a,v)/dE = ^dn{r, a)/dE 

a a 



(50) 
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Figure 7. Decomposition of the local density of states of a two terminal conductor into 
emissivites. 

Equation (50) represents a decomposition of the local density of states into 
emissivities and inj activities. We mention that the emissivities and injec- 
tivities are not independent of one another. In fact, in the absence of a 
magnetic field they are equal, dn{a,r)/dE = dn{r,a)/dE. In the presence 
of a magnetic field the microreversibility of the scattering matrix implies 
that the emissivity into contact a in magnetic field B is equal to the injec- 
tivity of contact a if the magnetic field is reversed, 

dn{a, r; B) /dE = dn{r, a; -B) jdE . (51) 

While the local density of states at equilibrium is an even function of the 
magnetic field, i. e. dn{r; B)/dE = dn{r; —B)/dE, the injectivities and 
emissivities are in general not even functions of B. This has some pecu- 
liar physical consequences, as has been shown recently by a low-frequency 
measurement of capacitances in a quantum Hall system [45] . 

5.5. COMBINED EXTERNAL AND INTERNAL RESPONSE 

We need to find an expression of the current response generated by the elec- 
tric potential oscillations caused by the external potentials Sjiai^) exp(— ia;t) 
We are interested in the response to first order in w, and since the currents 
in Eqs. (39) and (43) are proportional to oj it is sufficient to know the quasi- 
static nonequilibrium state discussed in Sect. 3. We express the deviation 
of the potential away from the equilibrium potential with the help of the 
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characteristic potentials defined in Eq. (11). ^From Eq. (43) we find 

6li{uj) = iew^ / (fr {dn{a,r)/dE) up{r)dnp{u) . (52) 

Thus the induced potentials give rise to a conductance 

d^r{dn{a,r)/dE)uf3{r) , (53) 
which can be written in terms the Green's function and the injectivity, 

Gipico) = ie^LO I d^r [ d^r'{dn{a,v)/dE)g{r,Y'){dn{Y' ,P)/dE) . (54) 



Eq. (54) tells us that the internal response G'^^g{uj) is a consequence of the 
charge injected from contact /? which generates a potential determined by 
the Green's function and that this potential in turn generates a current at 
a determined by the cmissivity into that contact. The total response is the 
sum of the external response (40) and the internal response (54), 

G^p{u)=Glp{u) + G],p{u:) . (55) 

We express it in the form 

Ga/3(a;) = G^piO) - iujE^p + 0(0;^) (56) 

and call EajB the (screened) emittance of the conductor. It is given by [5] 

E„,^e-'-^-e^l^rjd^r''^s(ry)^-^. (57) 

Before we apply this result in the next section to a few simple problems, 
we want to demonstrate that Gap is indeed current conserving, i. e. that 
the rows and columns of the emittance matrix add up to zero. Consider the 
first column. If we add En and £'21 the first terms in the emittance give 
the total charge dNn/dE + dN2i/dE injected from contact 1. In the second 
term the two emissivities add to give the local density of states. Now Eq. 
(22) is used. What remains is the integral over the entire volume of the 
injectivity which is just the total injected charge. Thus for a two terminal 
conductor in electrical isolation the emittance matrix satisfies E^ = En = 

—E12 = —E21 = E22- 



6. Examples 

6.1. EMITTANCE OF A METALLIC DIFFUSIVE CONDUCTOR 

Let us consider a mesoscopic metallic conductor connecting two reservoirs. 
In a metallic conductor charge accumulations are screened over a Thomas- 
Fermi screening length (apart from miniscule and more subtle Friedel-like 
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long-range effects [37]). If we assume in addition that the density varies 
not to rapidly then the local potential is directly determined by the local 
density. The local potential 5U{r) generated by an injected charge (5nj„(r) is 
determined by {dn{v) / dE) eSU {r) = 5nin{r). This corresponds to a Green's 
function which is a delta function in space and with a weight inversely 
proportional to the local density of states g(r, r') = (dn(r) /dE)~^S{r — r'). 
Using this in Eq. (57) gives an emittance in terms of densities only [19], 



There are no electric field lines outside the conductor. 

The wire with cross-section A ranges from x = —L/2 to x = L/2. The 
mean distance between the impurities is I. The reflection and the trans- 
mission probability per channel are R = 1 — J/L and T = l/L, respec- 
tively. The partial densities of states of reflected carriers are dNu/dE = 
dN22/dE = (1/2)(1 - l/L)(dN/dE) where dN/dE is the total density of 
states in the volume J7 = AL. The partial densities of states of trans- 
mitted carriers arc dN 12 /dE = dN2i/dE = {l/2){l/L){dN/dE). The dif- 
fusion equation for the diffusive metallic conductor implies the ensemble 
averaged and over the cross section averaged injectivities dn{x,l)/dE = 



{l/2L){dN/dE){l-2x/L) 33iddn{x,2)/dE = (1/2L) (dN/dE) (l + 2x/L). 



In the absence of a magnetic field, the cmissivitics are given by the same 
expressions. The linear dependence of the injectivities gives an ensemble 
averaged potential which drops also linearly. With the help of these expres- 
sions, we find for the emittance [13] 



Already for metallic wires longer then 3/ the emittance is a positive quan- 
tity. Since typically l/L is much smaller than 1/3 the metallic diffusive wire 
at low frequencies responds like a capacitor with an ensemble averaged ca- 
pacitance given by Eq. (3). 

6.2. EMITTANCE OF A PERFECT BALLISTIC WIRE 

Consider a ballistic wire of length L with N quantum channels. We apply 
again Thomas-Fermi screening. However, this is not very well justified and 
permits to obtain an estimate only. We also ignore the variation of the 
potentials near the contacts, which in a more realistic treatment might well 
give us a capacitive contribution. Each quantum channel contributes with 
a density of states 2L/hvn, were Vn is the channel velocity evaluated at the 
Fermi energy. The total density of states per spin is dN/dE = J2n "^L/hvn- 
In the absence of backscattering it holds dNn/dE = dN22/dE = and 




(58) 



E^ = {l/Q-l/2L){dN/dE) . 



(59) 
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dN2i/dE = dNi2/dE = (1/2) dN/dE. The injectivities dn{x,l)/dE = 
{l/2L){dN/dE) and dn{x,2)/dE = {l/2L){dN/dE) are independent of 
the space coordinate. This corresponds to a potential in the balhstic wire 
which is constant and midway between the electrochemical potentials at 
the contacts. Hence, using Eqs. (57) we find that an ideal perfect wire 
has a negative emittance given by Eq. (4). A ballistic wire responds like a 
conductor which classically is represented by a resistance and an inductance 
in series, and where the emittance can be viewed as a kinetic inductance. 



6.3. EMITTANCE OF A RESONANT DOUBLE BARRIER 

As a third example we consider a resonant double barrier [19]. The scat- 
tering matrix is s^g = ((Jq/3 — ir^r^/A) exp(i(5Q + iSjg), where Sa/s is the 
Kronecker symbol and the 6a are phases whose energy dependence can 
be neglected compared to the rapid variation of the resonant denominator 
A = E — Ej. — e5U — iT/2 with F = Fi +F2. Here, Er is the resonant energy 
at equilibrium and e6U is the deviation of the electrostatic potential away 
from its equilibrium value at the site of the long lived state in the presence 
of transport. With the help of the injectances dN^/dE = (rQ,/27r|Ap) and 
the total density of states dN/dE = (F/27r|A|^) the partial densities of 
state can be expressed in the following manner: For a = /? we find 

dE 2 dE 2^ dE dE ' ^ ^ 

where the plus and the minus sign correspond to « = 1 and a = 2, re- 
spectively. Here, i? = 1 — T is the refiection probability. For a 7^ /3, one 
has 

(61) 



dNgp _TdN 



dE 2 dE 

The unscreened injectances (emittances) are found by integrating the in- 
jectivities (emissivities) over the volume of the localized state, i. e. over the 
well. They are given by 

dNia/dE + dN2a/dE = dN^/dE . (62) 

In a Thomas-Fermi approach the characteristic potentials Ua in the well 
are determined by {dN/dE)ua = dNa/dE. This gives Ua = Tq/F. Using 
Eq. (58) gives an emittance [19] 

2 {dN,/dE) {dN2/dE) ( rV2^^\ 
^'^ = dNjdE [ |AP J 

For a symmetric resonant tunneling barrier Eq. (63) simplifies and is given 
by Eq. (5) . At resonance the emittance is negative reflecting kinetic (induc- 
tive) behavior, it is zero at half-width of the resonance, and it is positive 
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(capacitive) if the Fermi level is more than a half width above (below) the 
resonant energy. Clearly, Thomas-Fermi screening is not very realistic for 
such a conductor. Moreover, the quantization of charge in the well might 
play a decisive role. Nevertheless, these considerations indicate the charac- 
ter of the results that a more realistic treatment might yield and hopefully 
stimulate work in that direction. 



7. Summary 

We have developed a self-consistent discussion of mesoscopic electrical con- 
duction. The determination of the electrical potential in the presence of a 
dc-current, although unimportant for the discussion of the dc-conductances 
itself, permits to discuss the conditions under which the dc-conductance for- 
mulae are valid, it permits to calculate the first nonlinear corrections of the 
purely linear response, and permits to find the ac-conductances to first or- 
der in frequency. We applied the results to some simple examples. 

We have emphasized the case of two terminal conductors in electric iso- 
lation. The theory permits, however, also to discuss the effect of nearby 
capacitors and gates and in fact provides a mesoscopic description of ca- 
pacitances [5, 46]. Some implications of a mesoscopic theory of capacitance, 
like Aharonov-Bohm oscillations in capacitance coefficients, the gate volt- 
age dependence of persistent currents have already been the subject of 
recent works [47, 48]. 

The self-consistent nature of electrical transport is a consequence of 
the long-range Coulomb interaction of carriers. A self-consistent descrip- 
tion must, therefore, tackle an interacting many-particle problem. Conse- 
quently, the results of such a theory depend somewhat on the sophistication 
that is used to treat the many-particle problem. Here we have used a sim- 
ple Hartree approach. Since density-functional theory is nothing but an 
improved Hartree theory it gives results which look formally very simi- 
lar to the results presented here [49]. A stronger modification of the results 
discussed here can be expected in situations where one must take the quan- 
tization of charge into account [48, 50]. 

The theory presented here, demonstrates that interesting results can be 
obtained by investigating nonlinearities and ac-conductances. The theory 
demonstrates that it is necessary to treat nonlinearities and the ac-response 
self-consistently to conserve both charge and current. We are confident 
that experiments will eventually demonstrate the close connection between 
electrostatic questions and nonlinearities and ac-response. 
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